Precision measurement of the metastable 3 P 2 lifetime of neon 
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The lifetime of the metastable 3 P2 state of neon has been determined to 14.70(13) s (decay rate 
0.06801(62) s" 1 ) by measuring the decay in fluorescence of an ensemble of 20 Ne atoms trapped in 
a magneto-optical trap (MOT). Due to low background gas pressure (p < 5 x 10 _11 mbar) and low 
relative excitation to the 3 D3 state (0.5 % excitation probability) operation only small corrections 
have to be included in the lifetime extrapolation. Together with a careful analysis of residual loss 
mechanisms in the MOT a lifetime determination to high precision is achieved. 

PACS numbers: 32.70.Cs, 32.70.Fw, 32.80.Pj 



In contrast to its importance for various active fields 
of research covering such a wide range as atomic physics, 
quantum optics, and nuclear physics, there still has been 
no measurement of the natural lifetime of the metastable 
3 P2 state |jj of neon with sufficient precision. A selection 
of research activities profiting from an improved measure- 
ment include the quest for Bose-Einstein condensation 
(BEC) of metastable neon, advanced atomic structure 
calculations, the investigation of ultracold collisions, and 
even precision tests of the electroweak theory being cur- 
rently persued by investigating the nuclear decay of an 
optically trapped sample of 19 Ne in the 3 P2 state 

Ongoing research directed towards BEC of metastable 
neon atoms in the 3 P2 state [|| Q clearly will benefit 
from an accurate knowledge of the state's lifetime. This 
includes the optimization of the production process as 
well as the study of collision processes such as Penning- 
ionization and elastic s-wave scattering. Exciting new 
hysics complementing the work on metastable helium 
can be expected for a 3 P2 neon condensate. We would 
like to point out the study of higher-order correlations 
P| and the intriguing possibility of a modification of the 
^2 decay rate due to the high phase-space density or the 
phase coherence in a BEC ||. 

Significant advance in atomic structure calculations 
has already been and will be further triggered Q by 
the work presented here. To date only a preliminary ex- 
perimental value of 22 s for the 3 P2 lifetime of neon ex- 
ists [Tl| ] with no detailed investigations having been per- 
formed p2] |. A precise determination of the lifetime will 
put a more stringent test on theory. For most rare-gas 
atoms, the 3 P2 lifetime depends sensitively on electron 
correlations and relativistic corrections. The latter have 
only a minor effect for the case of neon, making neon 
specifically interesting for a critcal test of electron corre- 
lations. In addition, our measurement will close the chain 
of 3 P2 lifetime measurements for rare-gas atoms be yon d 
helium performed in recent years for Ar (38^ s [Q~3[ ) , 



Kr (39l4 s @ and 28.3±i;fs (m)), and Xe (42.9j^s 



and 42.41^3 s |Tj|), respectively. This will deliver the 



input for a detailed systematic investigation of the Z- 
dependence of metastable lifetimes. 

An accurate experimental determination of the 3 P2 
lifetimes of rare-gas atoms other than helium has only 
been achievable by the ability to prepare cold atom sam- 
ples in a MOT Q [jj, |1- For a precision lifetime 
measurement, the coupling of the 3 P2 to the 3 D3 state by 
the MOT light has to be considered since it will modify 
the observed decay rate. To date, two different methods 
making use of a variable MOT-on/off duty cycle have 
been employed, either by directly measuring the rate of 
UV photons produced by the decay during a MOT-off pe- 
riod 114, ll| or by extrapolating the observed decay rate 



for varying MOT-on/off duty ratios to vanishing MOT- 
on periods ]L3|. 

In our approach, we record the decay of fluorescence 
in a MOT (Fig. |l|) for different steady-state values of the 
population 7r 3 of the 3 D 3 state (0.005 < tt 3 < 0.47). This 
allows the extrapolation of the observed decay rate to 
7r 3 = thus giving the 3 P2 lifetime. Extreme care has 
been taken to achieve optimized experimental conditions 
such as low background gas pressure (p < 5 x 10~ n mbar) 
and MOT operation for such low values as 7r 3 = 0.005. 
This minimizes the amount of required corrections dra- 



in 10000 




60 

Time [s] 



120 



'Author to whom correspondence should be addressed. 



FIG. 1: Decay of the fluorescence for a 20 Ne ensemble in a 
MOT for the given excitations to the 3 D3 state of 713 = 0.005 
(A/r = -5.3, I/I = 1.75, p = 3.4 x 10 _11 mbar) and tt 3 = 
0.47 (A/r = -0.16, I/I = 24, p = 4.7 x 10 _11 mbar). Both 
sets of data consist of five curves each, demonstrating the high 
reproducibility of our decay measurements. 
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matically and the accuracy to which the 3 P2 lifetime of 
Ne can be determined is improved considerably. 

With no observable changes in shape and temperature 
of the atom sample, the temporal evolution of the ob- 
served fluorescence is proportional to the temporal evo- 
lution of the number of trapped atoms N(t) . The decay 
can be described by the differential equation 

^N{t) = -aN(t) - p Jn 2 (r,t) d 3 r + 0(n 3 ). (1) 

The experimental fluorescence data is fitted to a solution 
of this equation. No absolute calibration of the atom 
number is required for the determination of the decay 
constant a. In Fig. [l], non-exponential two-body decay 
due to intra-trap collisions can be seen during the first 
30 seconds. Two-body losses depend on the integral over 
the local number density n(r, t) and are described by the 
parameter /3, the absolute determination of which would 
require an absolute calibration of n(r, t), which is not the 
scope of this work. An uncalibrated value of (3, however, 
is determined in our fitting procedure. For all data runs 
no contributions of order n 3 or higher could be observed. 

At low atom numbers the decay is dominated by one- 
body losses with rate 

a = — (1 - tt 3 ) + — tt 3 +-/p + MC. (2) 

The parameter 7r 3 gives the population in the 3 D 3 -state. 
Thus (1 — 7r 3 )7V atoms decay from the 3 P2-state with life- 
time T2- The combined rate for all decay channels from 
the 3 D 3 -state other than the transition to the 3 P2-state 
(see Fig. |J) is given by Tg -1 . Background collisions con- 
tribute with a rate of jp ||l6f with p being the pressure 
in the vacuum chamber and 7 being a constant depend- 
ing on the atomic and molecular species present. The 
parameter A4C summarizes possible MOT losses. The 
main objective of this work is to determine the lifetime ti 
from the observed decay rates a by extrapolating to van- 
ishing excitation (7r 3 = 0), background pressure {p = 0), 
and MOT losses (M.C = 0). Achieving optimized start- 
ing conditions (small 7r 3 and p) and minimizing poten- 
tial contributions of AiC by a systematic analysis of the 
MOT characteristics have been essential. 
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FIG. 2: Detail of the level scheme of neon. Atomic levels are 
labeled according to LS- and Racah notation. Transitions are 
labeled by their multipolar character. 



Following fnj] , the population 7r 3 is determined by: 

,3 = 1 c -dlh y . ( 3 ) 

2 1 + C ■ I/Io + (2A/r) 2 

with the total intensity / of the MOT light fields, the 
saturation intensity Iq — irhcT/SX 3 — 4.08mW/cm 2 
(linewidth T = 2n x 8.18 MHz, A = 640.4 nm) and de- 
tuning A. The phenomenological parameter C accounts 
for the deviation of the effective total intensity in a three- 
dimensional MOT from a one-dimensional two-level sys- 
tem We adopt the notion of [|| and take C = 0.7 
with an uncertainty of 0.3. The intensity / is measured 
with an uncertainty of 10 %, the detuning A with an 
absolute uncertainty of 1 MHz. Therefore, for all data 
presented here the uncertainty in C is the dominant con- 
tribution to the uncertainty in 7r 3 . 

The experimental setup is shown in Fig. |^. A col- 
limated beam of metastable atoms is decelerated in 
a Zeeman-slower with the slowing laser beam passing 
through the center of the MOT. The MOT light field 
consists of six individual beams (diameter 22 mm) with 
spatially filtered beam profiles. All light fields are de- 
rived from a dye laser which is frequency stabilized to a 
neon RF-discharge. 

For each run, the MOT is loaded for 400 ms reaching 
a number of up to 4 x 10 s atoms, a maximum collision- 
limited density on the order of lxlO 9 cm -3 and a typical 
temperature below 1 mK. Then, the loading process is 
terminated by blocking the atomic beam and the slowing 
laser beam by mechanical shutters, and by switching off 
the magnetic field of the Zecman slower. In the following 
measurement period of two minutes, the fluorescence of 
the trapped atoms is recorded by a CCD camera cap- 
turing five image frames per second. Fluorescence de- 
cay curves are obtained by integrating the counts of each 
frame over the region where atoms are present. Great 
care has been taken to eliminate the influence of residual 
stray light. The decay curves span up to four orders of 
magnitude in signal. Fig. [j] shows two sets of data taken 
for different values of 7r 3 . For 7r 3 = 0.005 a decay rate 
of a = 0.06929(48) s _1 and a respective decay time of 
or 1 = 14.43 s is observed. 

In the following sections we describe our procedure to 
extrapolate the observed decay rate a to the 3 P 2 de- 
cay rate t^ 1 ■ For the extrapolation of a to vanishing 
background pressure a systematic investigation of the 
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FIG. 3: Experimental setup. 
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MOT decay as a function of pressure has been carried out 
(Fig. §(a)). The pressure is given as the nominal reading 
of the vacuum gauge (Balzers IKR070). The specific com- 
position of the background gas is not known. However, 
we could observe a linear dependence of a on pressure 
over a range of almost two orders of magnitude, indi- 
cating a pressure-independent composition of the back- 
ground gas. From a linear fit for data taken for < 0.1, 
we can determine the slope 7 = 5.2(1) x 10 7 mbar~ 1 s^ 1 . 
As a special feature of our setup, we can directly deter- 
mine a possible offset by simultaneously measuring the 
fluorescence and the rate of ions produced by collisions 
using a multichannel plate (MCP) ion detector. Under 
typical operating conditions we detect ions of different 
origin: (a) Ne + -ions originating from intra-trap collisions 
(rate oc Jn 2 (r, t) d 3 r). (b) Ions which originate from col- 
lisions of metastable neon atoms with residual gas con- 
stituents (rate oc N(t) x p). Thus, the ion count rate 
Rion(t) can be modelled by 

Rion(t) = CiN{t) +c 2 J n 2 (r, t) d 3 r, (4) 

where c± is assumed to be proportional to p. The val- 
ues of ci obtained for different pressures are shown in 
Fig. ^(b). From a linear fit we can determine the offset 
of the pressure gauge reading to (4 ±7) x 1CP 12 mbar also 
allowing for a vanishing offset. By combining the results 
for 7 and the offset, we can extrapolate the observed de- 
cay rates to p — 0. Due to the already low background 
pressure for our decay measurements only small correc- 
tions arise: At p = 5 x 10~ n mbar the background col- 
lision rate of jp — 1/(420 s) is more than a factor of 25 
smaller than the observed decay rate a. 

We also investigated the possibility of additional atom 
losses from the MOT (contribution M.C in Eq. (||)) at- 
tributed to a finite trap depth Uq or finite escape veloc- 
ity v e Q . We assume that a lossless MOT would result 
from an infinite trapping volume leading to infinite Uq 
and v e . In principle, this could be achieved for the ap- 
plied magnetic field gradient B' approaching zero if the 
MOT beams could be made infinitely large. In all practi- 



cal implementations the trap depth will remain finite due 
to the finite size of the laser beams. Thus, for a given 
beam geometry, a variation in B' should allow to investi- 
gate the MOT stability and to quantify the MOT losses 
A4C which include diffusive motion out of the trapping 
volume Jl9| leading to an exponentially decreasing loss 
rate with increasing trap depth and non-ionizing back- 
ground gas collisions |16| already being accounted for by 
the p — y extrapolation. The influence of the magnetic 
field gradient on the measured decay rate a is shown in 
Fig. U for two values of 773 = 0.016 and 0.41, respectively. 

In the most relevant case of small ir^ (Fig. |](a)), for 
sufficiently large B' > lOG/cm a variation of B' does 
affect the measured a not more than the uncertainties 
obtained from the fit. Thus we conclude that for our 
operating regime (shaded region in Fig. |(a)) the trap 
depth Uq and escape velocity v e are sufficiently large and 
that potential MOT loss contributions to the decay rate a 
are not significant. For large 7r3 (Fig. ||(b)), a variation of 
the MOT stability with B' is observed. Unfortunately, no 
quantitative model of the MOT at high saturation exists, 
the application of the loss models discussed above is not 
straightforward, and only a minimum and a maximum 
estimate of MC can be given: Since our decay data for 
large 1T3 have been taken at B' = 1.5G/cm which gives 
the best MOT stability, it is not clear whether MOT 
losses have a significant contribution on the determina- 
tion of a at all. On the other hand, a linear extrapola- 
tion of a to B' = should overestimate the corrections 
to be made 0, [y| giving MC = 0.005(2) s" 1 as an up- 
per limit. Losses for large ir^ only influence the 1^3 — > 
extrapolation which we take into account by incorporat- 
ing 0s _1 < MC < 0.005 s _1 into the uncertainties of this 
extrapolation. 

After this systematic study of the influence of back- 
ground collisions and MOT losses the determination of 
the 3 P2 lifetime T2 can be obtained by extrapolating the 
measured decay rates to vanishing population 773. Only 
data for a pressure below 5 x 10~ n mbar are considered 
and an extrapolation to p = is performed. The depen- 
dence of a p =o on n 3 is shown in Fig. |^ exhibiting a sim- 
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FIG. 4: (a) Decay rates a obtained for 713 < 0.1 after extrap- 
olation to tts — as a function of the nominal background 
pressure Pdisplay (b) Coefficients c\ proportionai to the num- 
ber of ionizing background collisions as a function of Pdisplay 
ailowing the determination of the vacuum gauge offset. 
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FIG. 5: Systematic investigation of MOT losses. Decay rates 
are piotted as a function of the magnetic field gradient for two 
different excitations (a) 7T3 = 0.016 and (b) 7T3 = 0.42. Shaded 
areas indicate the gradients of typical operating conditions. 
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FIG. 6: Decay rates a p= o as a function of the excitation 7T3 
allowing the determination of r 2 in the limit of 713 — > 0. 



TABLE I: Contributions to the r 2 = 14.70(13)s extrapolation. 



MOT decay at tt 3 = 0.005 


0.06929 s~ 


i 


±0.00048 s- 


i 


7T3 — > 0-extrapolation 


+0.00028 s" 


-l 


±0.00013 s~ 


-l 


p — > 0-extrapolation 


-0.00156 s- 


-l 


±0.00037 s" 


-l 


l/r 2 


0.06801s- 


l 


±0.00062 s" 


i 



pie linear behavior. Displayed uncertainties in a p=u are 
given by the combined uncertainties in the fit of a and 
uncertainties in the p — > extrapolation. To minimize 
uncertainties in the tt^ — ► extrapolation, only data sets 
for 7T3 < 0.05 and ir^ > 0.4 showing the smallest absolute 
uncertainties in tt^ have been included. The final value of 
the 3 P2 lifetime and decay rate are T2 — 14.70(13)s and 
t^ 1 = 0.06801(62) s, respectively . Table | gives a sum- 
mary of contributions and uncertainties for these values. 
Since stable MOT conditions are directly accessible for 
7r 3 = 0.005 and p < 5 x 10~ n mbar only unprecedented 
small corrections are required. The final value for only 
deviates from the uncorrected value of a -1 = 14.43 s by 



2 %. 

In addition to this result, extrapolating a p =o to 7T3 = 1, 
that is to a hypothetical total population transfer to the 
3 D3-state, is also possible. From this, we gain the first 
experimentally determined lower limit of T3 > 59 s for 
the combined rate t^ 1 of decay of the 3 D3-state via all 
channels except the direct decay to the 3 P2 state. How- 
ever, this value depends sensitively on the A4 /^-correction 
applied for large 7r3 and only the stated lower limit is ac- 
cessible. 

To conclude we would like to compare our result of 
T2 = 14.70(13) s to theoretical values obtained by dif- 
ferent atomic structure calculations. Prior to our work, 
published theoretical values were given in [^0) to 24.4 s, 
which can be rescaled to 20.0 s using more accurate input 
parameters, and in pi] to 22 s. Significant remaining dis- 
crepancies initiated a reexamination of atomic structure 
calculations for the whole set of 3 P2 lifetimes of rare-gas 
atoms. A preliminary result based on a multiconhgu- 
ration Dirac-Fock calculation gives a value of 18.9 s for 
neon |[o|. A high degree of sensitivity to the included 
electron correlations has been found. Theoretical val- 
ues are still expected to be larger than the experimental 
value, since the inclusion of additional electron correla- 
tions should further reduce the calculated lifetime. Re- 
fined calculations are forthcoming fTpf , Most recent ab- 
initio calculations using MCHF methods give a lifetime of 
16.9 s |f-2| . Considering the combined progress in exper- 
iment and theory, we anticipate that the determination 
of the 3 P2 lifetime of neon may continue to serve as an 
important test for precision atomic structure calculations 
in the future. 

We thank Y.-K. Kim for valuable discussions and are 
grateful for financial support by the DFG within the 
Schwerpunktprogramm SPP 1116. 
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